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Abstract. A nontrivial space-time structure of the hot system created at RHIC 
is the defining aspect of the physics of relativistic heavy ion collisions. Fem- 
toscopy through pion intensity interferometry provides direct access to the dy- 
namic geometric substructure of the freeze-out stage, and appears to confirm 
{p^ , source evolution via collective flow. Since flow is usually considered to be a bulk 

■ phenomenon, it is surpring to find a simple scaling in preliminary HBT radii 

l/"} \ from p+p, d+Au and Au+Au collisions. Investigating the light system data in 

detail, we discuss a new way to visualize the "fine structure" of 3D correlation 
^ , functions and the potential importance of long range correlations. 
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1. Introduction 



A central goal in the study of relativistic heavy ion collisions is to probe an equili- 
brated bulk system of partonic matter. Copious results obtained at the Relativistic 
Heavy Ion Collider (RHIC) indicate the early ideas of the creation of a gas of quasi- 
free quarks and gluons (the Quark Gluon Plasma or QGP) need to be replaced a 
non-hadronic "liquid" [ I] of strongly interacting color matter (the sQGP [|2])> in 
which the degrees of freedom are composite quasi-particles [[I]. Understanding the 
nature of this matter at the most extreme conditions is of fundamental interest. 

The strongly interacting nature of the sQGP is indicated by very strong col- 
lective flow [ |3] and energy loss of hard-scattered partons; these are indisputably 
the most trumpeted observations at RHIC. However, these exciting probes deal 
exclusively with momentum space, and thus fail to probe a crucial part of the story. 

As the apparatus, terminology, and techniques of modern heavy ion physics 
resemble more and more those of particle physics, it is worthwhile to remember that 
the distinguishing feature of our field is a focus on non-trivial space-time structure. 
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The defining aspect of the new state of matter hoped for at RHIC is a change in the 
space-time length scales characterizing the system and color correlations. Further, 
this new phase of matter can only be generated in collisions of "large" systems. 
Geometric features of the nuclear overlap in the initial state (especially at finite 
impact paramter) drive the further evolution of the system. In the intermediate 
state (during or after thermalization) , collective bulk behaviour develops, generating 
non-trivial correlations in momentum and coordinate space. Further, in the final 
state, anomolously large spatio-temporal scales may signal the the phase change 
from the color-deconfined to confined (hadronic) state. Clearly, space-time aspects 
of the system evolution are of fundamental importance to understanding the system 
at RHIC. In such a context, disappointing indeed are recent tendencies [ |2] to 
dismiss measurements (annoying in that they disagree with calculations) with strong 
connection to the space-time structure of heavy ion collisions. 

Two-particle correlations at low relative momentum (for historical reasons often 
called "HBT correlations" when performed with identical mesons) represent the 
most direct method to extract space-time scales from a femtoscopic dynamic source. 
Pion HBT has been used extensively to extract the size [0] E] and shape [E] of the 
freeze-out configuration at RHIC. In the next Section, we very briefly review the 
status of our understanding HBT in heavy ion collisions. This motivates the present 
study which compares small collision systems (p+p, d+Au, and peripheral Au+Au) 
to the largest ones (central Au+Au) . Preliminary results suggest a surprising scaling 
of length scales between these systems. In an effort to understand the signal for the 
small systems, we introduce a new representation of the correlation function which 
points to the importance of long-range and not yet understood correlations. 

2. Comforting femtoscopic systematics 

Heavy ion physics observables are always determined by an interplay between var- 
ious underlying processes. Thus, extracting meaningful information from any ob- 
servable (e.g. flow, spectra) requires systematic study of the observable over a large 
range of conditions in which the relative contribution between the competing physics 
changes. The history of heavy ion physics is littered with apparent understandings 
of underlying physics, based on agreement of theory with data under very limited 
conditions. But when the data is studied systematically (e.g. by varying collision 
energy), the theory is seen to fail. One hopes that we are not being similarly mis- 
led at RHIC, where several successful theories are considered valid only at these 
energies, and only for relatively central collisions, so that failure to reproduce data 
under different conditions is not considered a problem. 

Perhaps even more than other observables, femtoscopic measurements need to 
be systematically studied, since questions of "puzzles" in the space-time sector at 
RHIC have been raised. Here, we touch on some of the main femtoscopic systematics 
at that AGS, SPS, and RHIC. A full discussion may be found in [EJ. 

To confirm the most basic premise of interferometry, it is important and com- 
forting to note that transverse (to the beam) length scales increase monotonically 
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Fig. 1. An indication of the state Fig. 2. Pion HBT radii from central Au+Au col- 
of femtoscopy 15 years ago. Data lisions at RHIC, fitted [E3j with the Blast-wave 
from [E] compiled in [Ej. model [|TT] using flow and temperature parameters 

from momentum spectra. 

with the size of the colliding system. Fifteen years ago, this was almost all that 
was known about HBT radii and even that with very large error bars. Fig- 
ure H represents the status of systematic femtoscopy in heavy ion physics prior to 
initial runs at AGS and SPS. In the past decade-and-a-half, data and techniques 
have evolved tremendously, generating an equally tremendous range of systematic 
femtoscopic studies. 

The extracted transverse shape of the system at freeze-out also makes intu- 
itive sense. At AGS energies, for which elliptic flow (preferred expansion in the 
reaction plane) is weakest, the freezeout shape approximates that of the overlap 
region between the target and projectile. At RHIC, strong elliptic flow leads to a 
considerably rounder shape. 

Kinematic systematics encode more detailed and more important dynamically- 
generated geometric substructure. Rapidity dependent femtoscopy supports a freeze- 
out configuration which is approximately boost-invariant at midrapidity, underscor- 
ing the role of longitudinal flow. More interesting, however, are transverse collective 
flow components, since they must be generated in the evolution of the hot system 
and provide some information on thermalization. The strong decrease of transverse 
HBT radii (from pions and other correlated particles) with ttit is believed to result 
directly from space-momentum correlations induced by flow. 

HBT radius dependences on p? (c.f. Figure [2J are described by collective and 
thermal dynamics as are momentum-space observables in flow-dominated "Blast- 
wave" models [ II 1| . The picture behind such models (which approximate hydro- 
dynamical freeze-out scenarios) is illustrated in the left panel of Figure [21 There, 
arrows represent average boost velocities of local elements of the emitting system. 
Particles emitted with some p are preferentially emitted from one region (called 
the "homogeneity region" [E]) of the source; higher-pr particles are emitted from 
ever-smaller homogeneity regions. HBT radii reflect the size of this region. 

In principle, temperature gradients [ 13. or a non-flowing but cooling source 
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Fig. 3. Left: Global space-momentum correlations in a collectively-expanding 
system. The homogeneity region for particles with a given px is shaded. Right: 
Local space-momentum correlations correspond to a "super-" system composed of 
several subsystems, each identical to that in the left panel. The homogeneity region 
in this case is the sum of all the shaded areas- that is, essentially the entire "super- 
" system size. While the geometric configurations of the two scenarios are quite 
different, they would appear identical in purely momentum space. 

would also produce size-momentum correlations and HBT radii which drop with px- 
However, non-identical particle correlations [ 114) show clear evidence for position- 
momentum correlations (note that the average emission point depends on p in the 
left panel of Figure 0, which arise only in the flow-dominated scenario. 

3. Puzzles, puzzles 

Into the comforting story described above, we inject some unpleasant facts. Hydrodynamics- 
the simple theory assuming only thermalization and using an equation of state- 
does strikingly well in the description of the bulk momentum distribution of emit- 
ted particles. However, hydrodynamical calculations fail to reproduce HBT radii. 
Doubting the validity of hydrodynamics or taking seriously purely hadronic dynam- 
ical models at RHIC [ I15| is almost anathema. Thus, the first "puzzle" is the failure 
of "reasonable" models to reproduce the HBT radii. It has been suggested that this 
puzzle is associated with timescales; the source evolution and decoupling may occur 
more rapidly in Nature than in models. 

While hydrodynamics- even in the momentum sector- only works at RHIC 
(and, it is supposed, higher) collison energies, the second femtoscopic "puzzle" arises 
from the violation of expectations based on very general grounds. As the collision 
energy is raised through the point at which QGP is barely created, additional en- 
ergy should go into the liberation of new degrees of freedom. At this point, the 
pressure would drop and expansion would slow due to the the "softest point" of the 
Equation of State. This was expected to lead to a rapid rise in freezeout length- 
and time-scales, measurable via femtoscopy. Disapponting and puzzling indeed is 
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Fig. 4. Preliminary HBT radii from Au+Au 
and d+Au collisions with ^/snn of various cen- 
tralities have been divided by those from p+p col- 
lisions at the same energy. The ratio is plotted as 
a function of pr- 



Fig. 5. Cartesian ("out-side- 
long") projections of prelim- 
inary correlation functions 
from peripheral d+Au colli- 
sions measured by STAR. 



the observation of almost constant HBT radii as the energy scale is varied over two 
orders of magnitude 

When "puzzles" accumulate, it is well to revisit those aspects of the issue which 
are assumed to be well-understood. As mentioned above, the py-dependence of 
HBT radii is thought to encode important dynamical information; the homogeneity 
lengths decrease with increasing pt due to collective transverse flow. Bulk flow 
requires substantial reinteraction of system components- specifically, a mean free 
path much smaller than the system size- in order to develop. Indeed, the failure 
of hydrodynamics to successfully describe data from lower energies or peripheral 
collisions is attributed to insufficiently large matter density and/or system size. It 
might seem, then, that the collective flow signal in HBT radii from central Au+Au 
collisions, would be smaller in peripheral collisions and smaller still in p+p collisions. 

In prior measurements [ 9 , HBT radii from elementary collisions have been ob- 
served to decrease with pt, but this behaviour is generally attributed to resonances 
or string dynamics. STAR has also reported the decrease with pt in preliminary 
HBT results in p+p collisions [US]. Given the rather different underlying physics 
(bulk flow versus resonances or strings) behind the space-momentum correlations, 
the observation [^j of a trivial scaling R p+p (pt) ~ R Au+Au (pt) is striking. Fig- 
ure 01 shows the ratio of HBT radii from Au+Au and d+Au collisions of various 
centrality, divided by the HBT radii extracted from p+p collisions. Although this 
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scaling may be simply a benign coincidence, for now we call this the third femto- 
scopic "puzzle." 



4. Viewing and diagnosing fine structure in 3D correlation 
functions 

As is commonly done, our HBT radii are extracted via fits to a functional form 
which assumes a Gaussian distribution in space for the homogeneity region. Not 
surprisingly, the distribution is never found to be purely Gaussian in practice, and 
the consequences of non-Gaussianness (usually within a model-dependent context) 
have been discussed extensively in the literature. One of the few quantitative exper- 
imental evaluations of systematic errors associated with non-Gaussian components 
has recently been published [ 5 . There, we conclude that while such components 
are significant, the Gaussian HBT radii do characterize the correlations, even when 
a more general functional form is used. Here, we discuss a different aspect of the 
measured correlations for the lightest systems. 



4-1. Problems in (and with) Cartesian projections 

Correlation functions are measured and fit in three-dimensional relative momentum 
components q — q s — qi ("out-side-long") space [E| using a functional form [E] 

nt \ A {qo;is,qi) r, 1 J, , - q 2 R 2 - q 2 R 2 ~ q fRf\ 

C{q ,q s ,qi) = — = N ■ (1 - A) + \K cou i{q inv ) 1 + e q ° <• 1. a ii 1 

(1) 

where K cou i(qi nv ) is the squared relative wavefunction integrated over a Gaussian 
source of finite size, and N is a normalization factor independent of q. 

Visual evaluations of the quality of the data and fit are typically done via 1- 
dimensional projections of the data and the fit onto qi (i = o, s, I): 

, , £ 3 -,fc^(go,ft,gQ , Y. j ,kT{q ,q a ,q l )-B(q ,q a ,qi) 

Here, T(q a ,q s ,qi) represents the right-hand side of Equation^ Sums in the pro- 
jection for each component i are done over bins with qj and qk (k,j 7^ i) in some 
small range- typically on the scale of the fcmtoscopic correlation enhancement. 
Equation [3 defines what we call "Cartesian projections" in the out-side-long space. 

Figure[5]shows such projections for preliminary correlation functions for periph- 
eral d + Au collisions. Clearly, the fit and data projections in q a are significantly 
different in the femtoscopic (low-g) region. To first order, this suggests that the 
extracted R is too low by ~ 30%. However, this is not the full story, as we do 
not see the full 3D structure of the correlation functions in projections from Equa- 
tion |2J it may be that nontrivial structure in unplotted g-bins cause the fit (which 
globally minimizes some figure of merit such as negative-log-likelihood) to make 
compromises which show up in this particular projection. Similarly, non-experts 
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sometimes scratch their heads when viewing 1-dimensional projections: there is 
often no apparent reason why the fit should not do much beter than it does. 

In Figure however, we see clear problematic structure even in the Carte- 
sian projections for d+Au data. In particular, the correlation function in the 
large-q ("normalization") region has different values in the different projections. 
Pions emitted from p+p collisions have a very large Bohr radius compared to the 
source size, so structure at very large q cannot be of femtoscopic origin as quantum 
statistics-induced correlations die out, and Coulomb final-state interactions will be 
overwhelmingly s-wave. It could in principle arise due to global or other correla- 
tive effects such as momentum conservation or jets, though our preliminary studies 
suggest that this is not the case here. We do not present these studies or further 
discuss possible origins of the structure, but only focus on its characterization here. 
This preliminary result in Figure using low-multiplicity data from STAR repre- 
sents the worst-case dataset, and so is a good dataset on which to focus efforts to 
study nontrivial structure in the correlation functions and the extent to which they 
influence the extracted femtoscopic scales. 

4-2. Representation via spherical decomposition 

In principle, one could visualize the full structure of the 3-dimensional correlation 
function via a series of Cartesian projections in qi over different ranges in q^k- This 
would, however, constitute a large number of figures, and relevant patterns which 
cut across projections might not catch the eye. Given the symmetries of the corre- 
lation in g-space, a more natural representation is through spherical decomposition. 
Here, the spherical coordinates 9, 4>, and Q relate to the Cartesian ones as 

q = Q sin 9 cos 0, q s = Q sin 9 sin <f>, Qi — Q cos 9. (3) 

The correlation function may be binned 1 in Q, cos 9, and <fi and represented via 
its spherical harmonic coefficients, which will depend on Q: 

AxAQ)= J2 C{Q 1 9. ll <j )l )Y x ^(9. ll( j )l )F x ^(9 l ,A cos g,A ri> ), (4) 

bins i 

where F\ jfi represents a numerical factor correcting for finite bin sizes A cos g and 
A^; it turns out not to depend on c/>j. 

Symmetry constrains the number of relevant components. For femtoscopic anal- 
yses of identical particles at midrapidity which integrate over reaction-plane orien- 
tation (i.e. almost all analyses to date), only A\ tli 's with even values of A and 
\x do not vanish. Further, it is natural to expect that the statistical relevance of 
high-A components is diminished. Thus, by glancing at only a few one- dimensional 
plots, one views the entire correlation structure in orthogonal components. Once 
the practitioner is familiar with this representation, it should have diagnostic value. 

lr To avoid artifacts, bin widths should be = j^- and A cos g = — . 
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Another advantage of the spherical representation is that they are less sensitive 
to detector acceptance effects than arc the Cartesian projections. While the three- 
dimensional correlation function is ideally independent of acceptance, the summing 
and weighting with bins in Equation will affect C' l (qi). Thus, while NA44 and 
NA49 should measure the same correlation function in (f-space, their f-dimcnsional 
projections will differ. Their A>, M 's, however, will be identical, assuming that both 
experiments have some nonvanishing acceptance in all (f bins. Related to this, the 
one-dimensional correlation function C(Qi nv ) depends heavily on acceptance, while 
Aq q(Q) does not. Finally, from a physics standpoint, it has recently been pointed 
out [ I18j that A\ t n components correspond one-to-one with anisotropics of the same 
order A, fi of the spatial homogeneity region. 

To begin calibrating ourselves in this new representation, we show in Figure El 
harmonic coefficient functions corresponding to a purely Gaussian homogeity region. 
(In this simple example, K cou i in Equation ^ has been set to unity.) Some simple 
intuitive trends are clear. For this case, ^4o,o(Q) depends only on R? Q + R 2 S + Rf 
and not on shape information. The sign and magnitude of ^2,0 reflects the relative 
magnitude of the longitudinal to the transverse radii, independent of R Q /R S . In 
contrast, the important comparison of R D and R s is reflected in the sign of ^2,2- 
Information encoded in higher-order components is more complicated, but clearly 
simple degeneracies and sign-changes appear there, as well. 

Generic features are also clear. We notice further that, with increasing A, 
the "bump" in A\ M (Q)becomes increasingly smaller in magnitude, as hoped (so 

that only a few A terms are relevant). Also, A\^{Q) ^— >° if A 7^ 0, and the 
"bump" moves to higher Q with increasing A, reminiscent of wavefunctions for 
central potentials. Finally, the calculated correlation function approaches a constant 

value (here, unity) for large Q, independent of direction (6, <p)\ thus, A\^ — >°° So t \. 

The spherical decomposition for the worst-case (low multiplicity) d+Au data 
is shown in Figure together with the decomposition of the fit. A\ ttl components 
for A > 4 were statistically insignificant. Clearly, the functional form in Equa- 
tion [T] ( now with the Coulomb term) does not describe the correlation anisotropy 
or strength. The overall size probed by ^4o,o is probably accurately extracted, but 
the anisotropic shape (e.g. R /R a ) may be underestimated in the fit because of 
the influence of large-Q correlations. The nonfemtoscopic A ^ components ap- 
pear roughly independent of Q, and might contribute as well into the femtoscopic 
(low-Q) region. At present, we are actively investigating the source of the large-g 
anisotropy in the data and how to account for it in the fit. 

5. Summary 

We have presented preliminary HBT radii for d+Au and p+p collisions and pointed 
out a surprising scaling effect when they are compared with those published for 
Au+Au collisions. A naif interpreting Figure 01 would conclude that the linear scale 
of a gold nucleus is about six times that of a proton; despite the tremendous detail 
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accummulatcd since then, this threatens to return femtoscopy to its status of 15 
years ago (c.f. Figure With arguments of flow and space- momentum corre- 
lations and expansion 2 , together with a dose of coincidence, a more sophisticated 
expert would claim that much more information is encoded in the HBT systematics. 
However, Occam's razor is a formidable adversary. It is urgent to study the "third 
femtoscopic puzzle" in detail. 

We are exploring a new representation of the multidimensional correlation func- 
tion which has several technical advantages and a simple connection to physics. This 
representation makes visibly manifest the full structure of the correlations, and has 
revealed, for the small collision systems, long-range nonfemtoscopic correlations of 
unknown origin. At present, it is unclear the extent to which these not- understood 
features are related to the "third puzzle." Model studies into their origin, and 
attempts at more general fits to the data, are underway. 
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Fig. 6. Spherical harmonic coefficients A\^(Q) corresponding to a purely Gaussian 
correlation function with various HBT radii. 




3D Gaussian fit to the data 
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Fig. 7. Preliminary spherical harmonic coefficients A\ tll (Q) from peripheral 
d+Au collisions, together with the harmonics from a Gaussian fit. 



